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REPORT  SUMMARY 


At  the  beginning  of  the  quarter,  diagnostics  production 
and  the  diagnostics  generator  system  was  subcontracted  to  Automation 
Technology,  Inc.  of  Champaign,  Illinois. 

The  software  effort  was  directed  primarily  towards  preparing 
existing  codes  for  the  Biorroughs  B65OO  computer,  providing  a  coherent 
body  of  documentation  for  users  and  for  reference  purposes ,  and 
providing  a  COCKROACH  (ILLIAC  IV  FORTRAN)  to  GLYPNIR  translator. 
COCKROACH  programs  are  now  being  compiled  satisfactorily  and  work 
continues  on  extending  the  existing  subset  of  FORTRAN. 

During  the  quarter  an  interface  to  connect  the  PDP-11  System 
to  an  Interface  Message  Processor  (IMP)  was  specified  and  designed. 
The  PDP-11  simxilator  on  the  B5500  was  completed. 

Bids  were  received  and  final  purchase  orders  sent  for  hard- 
ware items  for  the  ARPA  Network  port  facility  for  the  Center  for 
Advanced  Computation.   Delivery  of  the  items  is  expected  in  January, 
1971. 

Negotiations  were  begun  with  Digital  Equipment  Corporation 
concerning  a  financial  support  agreement  for  the  development  of  the 
ARPA  Network  port  facility. 

Negotiations  on  the  final  cost  for  the  ILLIAC  IV  computer 
have  been  concluded.   The  result  of  the  negotiations  are  that  Burroughs 
Corporation  has  submitted  an  anticipated  cost-to-complete  of  27.7 
million  dollars .   Amendment  VI  of  the  contract  between  the  University 
of  Illinois  and  Burroughs  will  reflect  that  the  cost  for  the  computer 
system  will  be  27  million  dollars  and  that  Burroughs  and  the  Federal 
Government  will  enter  into  a  cost  sharing  agreement  up  to  31  million 
dollars . 
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The  contracts  for  the  Laser  Mass  Memory  System  were  signed  in 
late  December.   Precision  Instrument  Company  of  Palo  Alto,  California 
has  "been  contracted  for  this  system  which  is  to  be  attached  to  the 
ILLIAC  IV  computer  system. 

Project  expenditures  and  commitments  through  December,  1970: 

Burroughs  Corporation       $26,296,116.00 
University  of  Illinois        6,482,715-92 
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1 .   HARDWARE 


1.1  Diagnostics 

At  the  beginning  of  this  quarter  all  responsibility  for 
diagnostics  production  and  the  diagnostics  generator  system  was  sub- 
contracted to  Automation  Technology,  Inc.  of  Champaign,  Illinois. 
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SOFTWARE 


2.1  Operating  System 

2.1.1  Operating  System  I 

Diiring  this  quarter.  Operating  System  I  (OSl)  was  brought  as 
far  to  completion  as  the  B5500  would  allow.   Most  of  the  modules  have 
heen  run  together  and  have  been  found  to  work.   In  the  latter  half  of 
the  quarter  much  of  this  work  has  been  reproduced  on  the  B65OO.  The 
possibility  of  simiilating  ILLIAC  IV  through  its  own  operating  system  is 
being  studied. 

2.2  Compilers  and  Translators 

2.2.1  GLYPNIR 

GLYPNIR,  Version  II,  was  released  and  has  been  running  well. 
Version  II  includes : 

I/O  statements  for  the  simulator; 
Vectors ; 

Changes  to  some  control  statements  (notably  the 
<IF  statement>) . 

Version  II  was  also  converted  to  run  on  the  B65OO.   As  of 
December  31,  1970  the  compiler  works  correctly  on  the  B65OO  although 
the  6500  version  will  not  be  supported  until  the  end  of  the  first 
quarter  of  1971. 

2.2.2  COCKROACH 

A  program  which  translates  the  first  version  of  a  subset  of 
ILLIAC  IV  COCKROACH  (FORTRAN)  into  GLYPNIR  was  written  using  the  TWST 
translator  writing  system. 
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Some  minor  modifications  to  the  language  specifications  have 
been  made  to  enhance  translatability .   Several  benchmarks  will  nov  be 
written  and  tested. 

2.3  Assembler 

Work  began  on  translating  the  simulator  and  assembler  to 
the  B65OO,  but  this  is  proceeding  relatively  slowly  because  of  the 
difficulty  in  accessing  a  B65OO. 

A  collection  of  mathematical  macros  has  been  made,  debugged 
and  documented.   These  macros  (such  as  inner  (32  bit)  to  full  i6h   bit) 
precision  conversion)  together  with  trigonometric  and  elementary 
mathematical  subroutine  (SIN,  COS,  SQRT,  etc.)  will  form  an  excellent 
basis  for  a  library  of  functions  available  to  programmers  in  ASK, 
GLYPNIR  and  COCKROACH;  all  are  written  in  ASK. 

2 .h     Interactive  Communications  and  Graphics 

2.U.1  Interactive  Coimiiuni  cat  ions 

During  the  reporting  period,  several  major  software  efforts 
for  the  PDP-11  system  were  started,  completed  or  continued. 

•   The  B5500  PDP-11  simulator  was  completed. 

A  high-level  language  compiler,  PEESPOL  (runs  on  B5500), 
was  completed  as  Version  I. 

Ninety  per  cent  of  the  coding  to  be  done  on  the  ARPA 
Network  Terminal  System  (ANTS)  has  been  done — 50% 
debugged  by  simulation. 

A  paper  describing  software  development  for  the  B65OO  to 
II/EP  interface  was  completed. 

Bids  were  received  and  final  purchase  orders  sent  for  the 
following  hardware  items  for  the  ARPA  network  port  facility: 


-5- 


1  ea.  PDP-11/20-CA  system 

3  ea.  ^k  words  memory 

1  ea.  real  time  clock 

1  ea.  high  speed  paper  tape  system 

1  ea.  KSR-35  operators  TTY 

1  ea.  card  reader 

1  ea.  132  col-umn  line  printer 

1  ea.  256k  word  disk 

2  ea.  DECtape  drives 

12  ea.  2UOO  baud  line  interfaces 

2  ea.  150  baud  line  interfaces 

k   ea.  103A  dataset  controls  and  interfaces 

6  ea.  general  purpose  interfaces. 

Delivery  of  these  items  is  to  begin  in  early  Januaiy,  1971. 

An  interface  to  connect  the  PDP-11  to  an  IMP  was  completed  and 
specified.   Backplane  wiring  will  be  performed  by  Digital  Equipment 
Corporation  (DEC);  DEC  modules  will  also  be  used. 

Bids  were  received  on: 

Alphanumeric  CRT's  —  Four  terminals  have  been  ordered 
from  Hazeltine.   The  terminals  will  operate  at 
2UOO  baud. 

A  graphics  storage  scope  terminal  —  One  terminal  with  a 
joy  stick  has  been  ordered  from  Computek. 

An  electrostatic  printer  —  A  Gould  hQOO   system  is  now 
on  order  from  Automation  Techology,  Inc. 

Negotiations  were  begun  toward  a  financial  support  agreement 
with  DEC  on  the  development  of  the  ARPA  network  port  facility. 

2.1+.2  Graphics 

Aside  from  the  purchase  of  a  storage  scope  terminal  system, 
work  continued  in  the  area  of  contoxor  presentation  algorithms  and  surface 
function  (y  =  f(x,y))  presentation  systems. 
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2.5  Debugging  Package 

The  debugging  package,  SPIDER,  entered  its  initial  implemen- 
tation phase.   While  SPIDER  is  being  implemented,  the  program  DEVIL/BUGS 
is  available  for  debugging  through  the  SSK  simulator  on  the  B5500  and 
B65OO.   University  personnel  using  this  program  have  been  very  pleased 
with  it . 

2.6  B55OO/B65OO  Operations 

Throughput  on  the  B5500  was  very  good  during  the  quarter, 
although  a  few  days  were  lost  because  of  air  conditioning  problems. 

The  Burroughs  B65OO  computer  is  expected  to  arrive  at  the 
University  by  the  end  of  January,  1971.   Ideally,  the  system  will  be 
operational  by  the  end  of  February.   It  is  expected  that  the  B5500  will 
be  maintained  for  one  month  following  the  installation  of  the  B65OO. 

2.7  Software  Documentation 

A  tour  to  some  potential  ILLIAC  IV  users  and  a  user  question- 
naire have  helped  us  gain  insight  into  the  service  and  performance  users 
will  expect  of  ILLIAC  IV  software. 

An  ILLIAC  IV  Software  Reference  Manual  [l]  has  been  set  up, 
and  provides  a  basis  for  a  systematic  system  of  documentation,  a 
considerable  amount  of  which  was  provided  during  the  quarter. 
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3.   APPLICATIONS 


3.1  Numerical  Analysis 

3.1.1  Matrix  Inversion  and  Solution  of  Linear  Algebraic  Equations 

This  quarter,  the  Householder  triangularization  subroutine  for 
matrices  less  than  6h   x  6h   was  completed  and  is  executing  correctly  on 
the  sim\iLator.   This  routine  has  been  standardized  according  to  sub- 
routine linkage  prescribed  in  the  software  reference  manual  [l].   The 
Householder  method  allows  the  user  to  perform  several  matrix  operations, 
including  finding  the  inverse  and  solving  a  system  of  simultaneous 
linear  equations.   An  iterative  improvement  of  these  operations  is  also 
allowed.   Documentation  for  the  user  application  is  being  completed. 

Applying  the  Householder  subroutine  to  several  ill-conditioned 
test  matrices  revealed  that  double  precision  arithmetic  is  required  for 
resiilts  of  high  accuracy,  especially  for  large  order  matrices.   Present- 
ly code  is  being  written  for  the  general  size  matrix  in  single  precision 
and  provisions  are  to  be  included  for  disk  access  of  the  matrices.   The 
precision  routines  for  addition,  subtraction,  multiplication,  and 
division  written  by  Toshio  Yasui  [2]  are  being  tested  for  usage  within 
the  Householder  subroutine  for  the  general  matrix.   The  main  effort  now 
is  the  internal  computation  and  "checking  out"  the  particular  storage 
scheme  to  be  used  for  high  order  matrices. 

In  addition  to  standardizing  Householder  subroutine  linkage, 
the  method  of  modification  for  matrix  inversion  was  updated  to  comply 
with  these  standards;  documentation  is  being  completed  for  this  method. 

3.1.2  Eigenvalue  Problems 


3.1.2.1  Jacobi's  Algorithm 


Jacobi's  algorithm  for  finding  eigenvalues  of  real  symmetric 
matrices  has  been  written,  debugged,  tested  successfully  on  the 
simiilator,  and  documented. 
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The  Jacob!  algorithm  has  been  modified  to  accept  real  normal 
matrices  and  return  complex  eigenvalues  and  complex  eigenvectors.   This 
modification  has  not  been  tested. 

3.1.2.2  Eberlein's  Algorithm 

This  algorithm  has  been  coded  and  debugged.   It  is  presently 
being  tested  for  different  matrices.   Birring  these  tests  further  inves- 
tigations for  convergence  and  overcoming  oscillation  problems  were 
conducted  (see  [6]).   The  problems  arising  have  been  solved;  the 
solution  is  as  follows : 

In  normalizing  a  real  matrix  the  main  concern  is  to  maximize 

T    T 
C       ,  the  off-diagonal  elements  of  the  matrix  C  =  A(x)A  -A(x)A. 

This  is  accomplished  by  applying  orthogonal  transformation  matrices  to 

C  (see  [6]).   With  these  transformation  matrices,  the  new  C 

C        attains  its  maximum  value,  C   t  ov  ~  -'-/^  h,  where 


2k-l,2k 


^  "  ^^^2k-l,2k  ■"  ^^2k-l,2k-l  "  ^2k,2k^  ' 


1/2 


(1) 


1  2 

To  maximize  C^,  _  ^,  ,  therefore,  is  to  maximize  C^,  ^  ^,  and 
2k-l,2k'  2k-l,2k  


^  2k-l,2k-l    2k,2k^ 


So 


far  the  attention  [6]  has  only  been  on  maximizing  C 


2k' 


In  order  to  maximize  C  -  C     ,  let  us  first  look 

at  an  example : 

Assuming  the  diagonal  elements  of  C  are  of  the  form 


f-h 


^1.2 


'0.5 


1,2 


0.5> 
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we  see  that  C   -  C,  > 


=  0  and  it  will  not  contribute  to  maximizing 


C'l  .  The  transformation  matrix  pertaining  to  this  case  will  be  an 
identity  matrix  and  the  element  C^i  will  not  change. 

In  exchanging  row  2  and  row  3  by  means  of  an  orthogonal 
transformation,  C  becomes 


'-h. 


0.5> 


3. 


0.5> 


and  C__  -  C„,  =2.5.   Care  must  also  be  taken  that  the  C. .  elements  in 
33    34  ij 

the  upper  triangle  of  C  are  close  to  max|C. .|. 

i^j  ^^ 
Tests  have  shown  that  taking  both  C 


2k-l,2k'  ^""^   ^^2k-l,2k-l 


C^.  ^,  )  into  account,  oscillations  are  overcome,  rate  of  convergence  is 
2k, 2k  '  '  ^ 

faster,  and  the  result  (matrix  normalization)  is  more  accurate. 


3.1.3  Roots  of  Polynomials 

Newton's  Method  for  evaluating  square  roots  and  solving  non- 
linear systems  of  equations  has  been  coded  in  ALGOL  and  debxigged.  The 
method  for  evaluating  square  roots  has  also  been  coded  in  ASK. 

Work  on  Newton's  Method  in  solving  general  polynomial  equa- 
tions has  begun.   Preliminary  steps  in  the  general  root  Finder  include 
the  use  of  the  reciprocal  polynomial  for  finding  roots  outside  the 
unit  circle.   Algorithm  applicability  has  not  yet  been  decided  for 
complex  and  multiple  roots. 

3.2  Linear  Prograjmning 

Code  development  and  testing  have  continued  fairly  smoothly 
this  quarter.   System  flowcharts  have  been  produced  for  the  modified 
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Hellerman  and  Rarick  reinversion  algorithm  for  INVERT.   ILLIAC  IV 
assembly  code  (ASK)  has  been  written  for  about  10%   of  the  subroutines 
for  this  INVERT  procedure.   Testing  of  these  codes  will  be  initiated 
in  the  next  quarter. 

During  this  quarter  the  manpower  assigned  to  LP  has  gradually 
been  reduced  to  about  one-half  of  the  previous  level.   This  has  been 
accomplished  by  reassignment  of  LP  personnel  to  other  CAC  projects  and 
to  increased  responsibilities  of  LP  personnel  in  other  areas.   This 
reduction  is  in  keeping  with  the  progress  of  the  LP  code  and  the 
schedule  for  further  development.   It  is  felt  that  the  ASK  programming 
and  systems  experience  of  LP  personnel  will  be  of  value  to  other 
applications  and  software  groups. 

3.3  Long  Codes 


In  an  effort  to  improve  on  the  results  obtained  with  the 
previously  developed  system  identification  algorithms  based  on  the 
stochastic  approximation  approach  of  R.  C.  K.  Lee,  the  emphasis  has 
been  shifted  to  correlation  techniques  based  on  the  work  of  Anderson, 
et  al.  [3],  and  Mehra  [h] . 

Two  identification  schemes  based  on  the  work  of  Anderson, 
et  al.,  have  been  implemented  and  are  being  tested.   One  algorithm, 
closely  following  the  approach  of  reference  [3],  requires  that  the 
observation  matrix,  which  relates  the  state  vectors  of  the  system  to 
the  observation  vectors,  must  be  an  identity  matrix.   This  scheme 

approximates  the  system  transition  matrix  directly,  thus  requiring  the 

2 
determination  of  n  values.   A  second  algorithm  adapts  this  approach 

to  the  formulation  used  in  the  previously  developed  stochastic  approxi- 
mation procedures  and  determines  an  estimate  for  the  n-vector,  i|f,  that 
represents  the  coefficients  of  the  characteristic  equation  of  the 
system  transition  matrix.   This  algorithm  should  give  results  compar- 
able to  the  first  algorithm  while  only  determining  n  values  instead  of 
n  values . 
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Work  with  the  correlation  approaches  indicates  that  it  may  be 
possible  to  extend  them  to  obtain  estimates  for  the  covariance  matrix 
of  the  observation  noise.   Since  knowledge  of  the  noise  covariance 
matrix  is  extremely  useful  for  identification  or  filtering,  the  possi- 
bility of  determining  observation  noise  covariance  is  appealing,  and 
work  is  continuing  in  this  direction. 

3.h     Signal  Processing 

A  GLYPNIR  program  has  been  written  which  computes  the  auto- 
correlation, c,  of  any  given  vector,  X.  : 

N-i 
c.  =  E  X.  .X.,    where  i  =  0,  1,  2,    ...,   N'  and  N'  <  N. 

This  program  is  considerably  slower  than  the  corresponding  ASK  program 
[5],  partly  because  GLYPNIR  cannot  access  registers  as  efficiently  or 
use  the  32-bit  mode. 

A  GLYPNIR  subroutine  is  being  written  which  will  compute  the 

fast  Fourier  transform  of  any  real  vector,  f.,  of  length  N  =  Gh .      This 

J 
subroutine  gives  the  real  and  imaginary  parts  as  well  as  the  amplitude 

of  the  Fourier  transform,  g.: 

N-1 
g  =  E  f  exp  (-2jtijk/N),     j  =  0,  1,  2,  . . . ,  N  -  1. 
^        k=0  ^ 

This  subroutine  will  be  extended  to  cover  any  length  which  is  a  power 
of  two  between  Gh   and  U096,  inclusive. 


3.5  Education 

3.5.1   ILLIAC  IV  Text 

The  initial  version  of  the  ILLIAC  IV  Textbook,  "An  Introductory 
Description  of  the  ILLIAC  IV  System"  will  consist  of  three  chapters — 
Backgro\and,  Hardware  Structure  and  Programming  Languages — and  will  be 
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issued  next  quarter.  Additional  chapters  will  be  issued  at  a  later 
time. 

3.5.2  ILLIAC  IV  Seminar 

A  one-day  seminar  on  ILLIAC  IV  was  held  on  November  3.      The 
outline  was  basically  the  same   as  the  outline  noted  in  the  previous 
QPR    [6]. 

3.5.3  Computer  Science-^91B 

The  graduate  course  "Architecture,  Applications  and  Languages 
for  a  Parallel  Computer"  will  be  offered  again  this  Spring  in  Room  237 
Digital  Computer  Laboratory  on  Monday,  Wednesday,  and  Friday  at  9  AM. 
The  outline  is,  again,  basically  the  same  as  the  one  presented  in  the 
last  QPR. 
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ADMINISTRATION 


k.l     Administration  and  Services 

if. 1.1  Contract  Status 

Negotiations  on  the  final  cost  for  the  ILLIAC  IV  Computer  have 
been  concluded.  Burroughs  has  submitted  an  anticipated  cost  to  complete 
of  27.7  million  dollars.   The  contract  between  the  University  of 
Illinois  and  Burroughs  Corporation  will  reflect  that  the  cost  for  the 
computer  system  will  be  $27  million  and  that  Burroughs  and  the  Federal 
Government  will  enter  into  a  cost  sharing  agreement  up  to  $31  million. 
The  schedule  for  the  cost  sharing  agreement  is  as  follows: 

Share  Ratio 

Amount Government /Burroughs 

Up  to  $27,000,000  100^  /  0% 

$27,000,000  to  $28,000,000  90%  /  10^ 

$28,000,000  to  $29,000,000  80^  /  20^ 

$29,000,000  to  $30,000,000  6o%  I  ko% 

$30,000,000  to  $31,000,000  10^  /  90^ 

The  scheduled  acceptance  of  the  ILLIAC  IV  is  J-une   1971.   The 
Center  for  Advanced  Computation  is  responsible  for  applications  work 
and  the  ILLIAC  IV  Project  is  responsible  for  providing  software  for  the 
ILLIAC  IV. 

During  December,  the  prime  contract  between  RADC  and  the 
University  of  Illinois  for  the  Laser  Mass  Memory  system  was  signed.   The 
subcontract  between  the  University  of  Illinois  and  Precision  Instrument 
Company  was  also  executed. 
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i|.1.2  Financial  Report 

Budgeted  expenditures- — Second  Quarter,  Fiscal  Year  1971: 


October 


November 


December 


Burroughs   $1,7^7,936.00 
University     215,302.00 


$552,3^0.00   $793,811.00 
215,302.00    215,302.00 


Total  for 
Quarter 

$3,09^,087.00 
6i+5, 906.00 


Actual  expenditures  and  commitments — Second  Quarter,  Fiscal  Year  1971 
(through  November,  1970): 


October 


November 


Burroughs     $l,10i+,  216.00*  $296,000.00 
University       22^1,663.50    195,011.56 


Total  through 
November  1970 

$1,400,216.00 
i+19,675.16 


Total  budgeted  and  actual  expenditures  to  date ; 


Burroughs  Corporation 
University  of  Illinois 


Total  Budgeted 

Expenditures 

through  December 

$26,296,116.00 
6,1+82,715.92 


Total  Actual 

Expenditures 

through  November 

$26,219,000.00 
6,266,182.95 


Includes  costs  of  ECS,  Texas  Instruments,  and  Fairchild  as  reflected 
in  recent  negotiations  with  Burroughs  Corporation. 
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